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ABSTRACT 

Context. One science goal of the upcoming Planck mission is to perform a full-sky cluster survey based on the Sunyaev-Zel'dovich 
(SZ) effect, which leads to the question of how such a survey would be affected by cosmological models with a different history 
of structure formation than ACDM. One class of these models are early dark energy (EDE) cosmologies, where the dark energy 
contribution does not vanish at early times. 

Aims. Since structures grow slower in the presence of EDE, one expects an increase in the number of galaxy clusters compared to 
ACDM at intermediate and high redshifts, which could explain the reported excess of the angular CMB power spectrum on cluster 
scales via an enhanced SZ contribution. We study the impact of EDE on Planck's expected cluster sample. 

Methods. To obtain realistic simulations, we constructed full-sky SZ maps for EDE and ACDM cosmologies, taking angular cluster 
correlation into account. Using these maps, we simulated Planck observations with and without Galactic foregrounds and fed the 
results into our filter pipeline based on the spherical multi-frequency matched filters. 

Results. For the case of EDE cosmologies, we clearly find an increase in the detected number of clusters compared to the fiducial 
ACDM case. This shows that the spherical multi-frequency matched filter is sensitive enough to find deviations from the ACDM 
sample, being caused by EDE. In addition we find an interesting effect of EDE on the completeness of the cluster sample, such that 
EDE helps to obtain cleaner samples. 

Key words. Galaxies: clusters: general - cosmic microwave background - Methods: numerical - Space vehicles: Planck 



1. Introduction 



Recent observations with ground based CMB experiments like 
BIMA, CBI, and ACBARlPawson e"tani2002t iFLeadhead et"an 
l2004t iReichardt et alj|2008l) report an excess in angular power 
at high multipoles { > 2000 with respect to the the- 
oretically expected CMB power spectrum. If explained by 
the contribution of the thermal Sunyaev-Zel'dovich (SZ) ef- 
fect jSunvaev & Zeldovichlll972h within the framework of the 
ACDM cosmology, a rather high normalisation of the mass fluc- 
tuations cTg ~ 1 would be required (Bond et al. 2005). Such a 
high value would contradi ct the erg = 0-796 ± 0.036 reported 
by the WMAP-5 results dPunklev et alj l2008i) . Furthermore, 
despite being controversial, the observed abundance of grav- 
itatio nal arcs (Bartelmann et al. 1998; Meneghetti et al. 2000, 
2003; Wambsganss et al...2004; Dalai et al. 2005; Hennawi et al. 



2007t iFedeli et al.1 l2008h and their presence in high-redshift 
galaxy cluste rs with z > 1 ('Hasinger et al."1998; Gladderj.eLalJ 
[2003; Thomp son et al. ! 2001; Zaritskv & Gonzalez 2003), gives 
rise to tension with the concordance model. 
One explanation of these contradictory observations is offered by 
models of early dark energy (EDE), in which the non-vanishing 
contribution of a dark-energy density at high redshifts alters 
structure formation. The effect of a few percent of EDE on non- 
linear stru cture formation has been analytically studied for in- 
stance by iBartelmann et al.l (l2006l) . who find that such a con- 
tribution may lead to ten times more galaxy clusters at redshift 
unity than ACDM. In this way the strong-lensing measurements 
(Fedeli & Bartelmann 2007; Fedeli et a l. 2008) and the ob served 
high-multipole CMB power spectrum dSadeh et al.ll2007l) could 
be explained. 



In addition, EDE models, belonging to the class of dynami- 
cal dark energy models, also offer an explanation for the present 
observ ed cosmic acceleration . Based on the concept of scalar 
fields ('Ratra & Peebles' '1988; 'Wetterich' '1988"; 'CaldwelleLal] 
1998; Liddle & Scherrer 1999; Zlatev & Steinhardt 1999) exert- 
ing negative pressure, these models are able to reproduce the 
late-time accelerated expansion of the Universe. The dynamical 
evolution of such a field, if present at all, is largely unconstrained 
by present observations of SN la, LSS and CMB and leaves room 
for a range of different classes of models (D pran e t al. 20051 
2007). M odels gi ving rise to attractor solutions iRatra & PeeblesI 
1988; , WetterichI IT9 8 8i) . where the evolution of the dark energy 
component follows the dominant component of the cosmic fluid, 
naturally suggest that dark energy might be non-negligible dur- 
ing long periods of the evolution of the Universe. 

Because a few percent of EDE significantly change the 
abundance of galaxy clusters with respect to the ACDM case at 
higher redshifts, we have reason to believe that such an excess 
could be detectable by the Planck mission. To investigate this 
we decided to simulate full-sky Planck observations of the SZ 
effect for models of EDE in order to quantify the impact on the 
expected cluster sample. 



Section 2 briefly reviews the concepts of EDE and is fol- 
lowed by a recapitulation of the SZ effect in Sect. 3. The con- 
struction of semi-analytic full-sky SZ maps for different cos- 
mologies is explained in Sect. 4. Section 5 outlines how realistic 
simulations of Planck observations are performed. The method 
of cluster extraction is discussed in Sect. 6, followed by an ex- 
ploitation of the obtained cluster sample for the considered cos- 
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Fig. 1. Properties of the considered cosmologies. The first panel on the left shows the redshift evolution of the equation-of-state 
parameter w{z) (left axis) and of the dark energy density Q.d^{z) (right axis). The middle panel shows the evolution of the growth 
function with scale factor a normalised to today (upper curves) and to early times (lower curves). The panel on the right shows the 
CMB power specti-a as computed with CMBEASY (,Doraa,2005i) . 



mological models in Sect. 
Sect. 8. 



7 and finished by the conclusions in Table 1. Key parameters of the used cosmological models. 



2. Early Dark Energy cosmologies 

2.1. General features 

Accelerated expansion in general can always be described by an 
equation-of-state parameter wde = Pde/Pde smaller than -1/3, 
either being constant as for example in the case of the cosmolog- 
ical constant wde = -1, or redshift depended wde - wde(z) as 
in the dynamical dark energy cosmologies. The special case of a 
non-vanishing dark energy density Qde(z) through long periods 
of the cosmic evolution can be achieved by allowing that wde(z) 
rises above zero. The left panel of Fig. [1] shows the redshift evo- 
lution of Q.de{z) and wue(z) for the different cosmological mod- 
els. For wde(z) the tracking behaviour is clearly visible for the 
EDE 3 and 4 cases and less pronounced for EDE 2. The nomen- 
clature will be introduced in Sect. 12.31 The energy densities for 
the EDE models do not vanish and are almost constant at z > 10, 
whereas for ACDM almost no contribution is left. A simple pa- 
rameterisation of EDE models has been introduced bv lWetterichI 
(|2004) and also by Doran & Robbers (2006) in which three cos- 
mological quantities are sufficient to fully determine the model; 
the present dark energy density Qde.o, the present equation-of- 
state parameter wde.o and the averaged value of the dark energy 
density during the era of structure formation. 



^DE,sf = -(Infleq) 



•J\na„ 



Q.UE{a)<i In a. 



(1) 



where cteq denotes the scale factor at matter-radiation equality. 
Current observational constraints from SN la, LSS and CMB 
data a llow for a Qde.si on the percent level dDoran et al. I I2OO5I 
l2007h . 



2.2. Impact on structure formation 

The study of structure formation in the p resence of dark en- 
ergy has been subje ct of ma ny studies (Wang & Steinhardt' 



19981: iDoran et al 



Mota & van de Bruck 



120011: Hinder & Jen kins 2003; 
I2004t iHorellou & Bered l2005l) . lead- 



Parameter 


ACDM 


EDE2 


EDE3 


EDE4 




0.265 


0.364 


0.284 


0.282 




0.735 


0.636 


0.716 


0.718 


i^DE.sf 




0.04 


0.033 


0.048 


h 


0.71 


0.62 


0.686 


0.684 




0.8 


0.78 


0.715 


0.655 


«s 


0.948 


0.99 


0.978 


0.976 


Wo 


-1 


-0.99 


-0.942 


0.934 



ing to the result that the growth of structures is quite sensitive to 



the presence of dark energy . I n the case of EDE t he non -linear 
part has been analysed by iBartelmann et alj (|2006'), who 
showed that at the high mass end a significant enhancement 
of the number of clusters at redshifts z ~ I can be expected. 
From the present point of view, the structures "decay" much 
slower towards the past than for the ACDM case. This could 
be a possible explanation fo r the un expectedly high lensing 
efficiency of distant clusters jFedeli & Bartelmann 2007) and 
the report ed excess power a t cluster scales in the CMB power 
spectrum dSadeh et al.ll2007h . 



2.3. The selected cosmological models 

For the present work we decided to study four cosmological 
models, three EDE cosmologies and the concordance ACDM 
model for comparison. To be consistent with p receding studies 
we an alysed one of the two EDE models of Bartelmann et al.l 
( l2006h namely the EDE 2 model. In addition we used two other 
EDE models (3, 4) kindly provided by Georg Robbers, both be- 
ing in agreement with all current bounds on EDE from SN la, 
LSS and CMB data dDoran & Robbersl2006l;lDoran et alj2007l) . 
where model EDE 4 resides at the limits, giving a bordering case. 
The key parameters for all models can be found in Table [T] For 
all figures we chose the same colour coding to distinguish the 
four different models, such that red illustrates the EDE 2 model, 
blue and cyan depict the EDE 3 and 4 cases and grey/black de- 
notes the fiducial ACDM case. 
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Cosmological model 

^m.'f^DE. ""S' ^^)^ 5d(^). ^vi^)' ^H' fga 
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spectrum P(/c,z) 



effective 
bias b^„{z) 



Limbers's 
\approximatiorT, 



clustered power 
spectrum Pyik.z) 
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I 
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of correlated clusters 
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using synalm 











I 



normalisation 



Poisson distribute 
clusters 
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cluster 
positions 



main routine: mass and redshift distribution 
► loop on AM 



loop on 



T 

■on re 

T 



Poisson realise 

N in (AM,Az) 



random realise 
M and z in (AM.Az) 



loop on clusters 



full-sky 
SZ map 




write HealPix map 
to file 



compute y-signal of cluster for given M and z 



compute total number 
of cluster electrons 




compute viral r„„ and 
core radius 
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electron density n^o 




uniform random 
position on sky 



compute pixel values 
scale to preserve Y 



store to HealPix map 



I 



if 1Cy* arcmin^ 
write to catalogue 



Fig. 2. Flowchart for the creation of full-sky Compton-y maps. The left column illustrates the process of obtaining the spatial cluster 
distribution on the sky, the centre one shows the routine for calculating the mass and redshift distributions, and the right column 
summarises the computation of the Compton-y profile for the individual clusters on the HEALPix sky map. Green boxes illustrate 
physical assumptions, blue ones computations, white ones stand for algorithmic structures and yellow for results of the routines. 



3. The thermal SZ effect 

The thermal SZ effect has in recent years become a popular tool 
for the search for clusters by current and future CMB experi- 
ments. The science case of the upcoming Planck mission also 
contains the task of delivering a full-sky cluster catalogue based 
on the SZ effect in order to study the cosmological implications. 
The effect is based on the inverse Compton scattering of CMB 
photons off the hot electron gas in the gravitational potential of 
the galaxy clusters. It is identifiable by its unique spectral sig- 
nature, such that at frequencies below v = 217 GHz the clus- 
ters are observed, with respect to the CMB, as shadows and 
above this frequency as emitting sources. The relative temper- 
ature change AT/T as a function of the dimensionless frequency 
X = hvlik^Tcus) can be calculated by 



AT / + 1 \ 

— -4, (2) 



where the amplitude, also known as the Compton-y parameter, 
is the line-of-sight integral of the electron pressure 

yiff) = ^ r nj.dl. (3) 

Here k^ denotes Boltzmann's constant, cti is the Thompson 
cross-section, mg is the electron mass, c denotes the speed of 
light, and and are the electron density and temperature, 
respectively. However, one is often more interested in the total 
Comptonisation Y, where y is integrated over the observed clus- 
ter surface 

y = r d^y{0) = £>a'(z)^ r dVne^e, (4) 
J m^c^ J 

which can also be rewritten as an integral over the cluster vol- 
ume, where Da denotes the angular diameter distance of the ob- 
served cluster at given redshift z- For the case of isothermahty 
the integrated Comptonisation is given by 

ksTs o-j 
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where A^e denotes the total number of thermal electrons within 
the cluster, being fully determined by the cluster mass Md via 



A'. = 



• 1+/h 
2mn 



/gasA/cl 



(6) 



for a known baryonic gas mass fraction /gas and hydrogen frac- 
tion /h, where Wp stands for the proton mass. For the first we 
assume throughout this work the cosmic value /gas = Qb/Qm - 
0.168 from WMAP-5 and for /h the generic value /h = 0.76. 



4. Creation of full-sky SZ maps 

In order to study the impact of EDE on the detectable cluster 
sample, it is necessar y to create full -sky maps of the thermal SZ 
effect. As outlined in lSchafer et al.l (l2006 a) one possibility is to 
use large scale dark matter only cosmological simulations in or- 
der to obtain the cluster positions and to additionally simulate in- 
dividual clusters including gas physics to get the Compton-y sig- 
nature, or t o model the individual c lusters by analytic j6-profiles, 
as done by iGeisbiisch et all (l2005h . Afterwards, everything can 
be projected on the sphere to obtain a full-sky Compton-y map. 
Since for EDE cosmologies these numerical simulations are up 
to now not available and in order to be more flexible, we decided 
to take a semi-analytic approach that is outlined in the following 
subsections and summarised in the flowchart in Fig. |2] For an 
overview on halo modelling in gener al the reader is referred to 
the review bv lCoorav & ShethI (l2002h . 



4.1. Collapse parameters for EDE models 

As mentioned in Sect. |2] the history of structure formation is 
quite different from the ACDM case. Basically two quantities 
are affected, the linear density contrast at collapse time 6c, and 
the growth factor For the calculations in the following subsec- 
tions these are crucial quantities, such that we briefly summarise 
their computation and properties. The dis cussion of the spherical 
collapse with EDE bv lBartelmann et all (|2006|) led to the result 
that the critical density contrast is given by 



lim 



£>+(xc) 



[A(x) - 1] 



(7) 



where x - a lathis the cosmological scale factor normalised to 
unity at turn-around, 0+ is the growth factor and A(jic) is the 
non-linear overdensity. The evolution of the growth factor for 
the EDE models is obtained by numerically solving the differen- 
tial equation 



+ Ifl E) ^ 2fl5£2 



£>+ = 0, 



(8) 



where E{a) is the expansion function of the Friedmann equation 
= hIe^, which is for the case of a time-dependend equation- 
of-state parameter w(a) given by 



E{d) - 1 0„fl + QoE exp 



do' 



-3 I 

h a 



1/2 



(9) 



Therefore all cosmological quantities depending on the expan- 
sion function are modified for EDE cosmologies. 
The virial overdensity Ay is onl y weakly affected b y the pres- 
ence of EDE as discussed in Bar telmann et a n (l2006h . where the 
way of calculating Ay is explained in more detail. 



4.2. The mass and redshift distribution of the cluster sample 

The computation of the mass and redshift distribution is the core 
part of the map construction and is summarised in the second 
column of Fig. |2] The basic information needed is the number of 
galaxy clusters per mass and redshift bin, that can be estimated 
using the comoving mass function 



dA^ dV dn 

-^{M,z) = AQ-^(z)-^(M,z), 
dMdz dzdO. dM 



(10) 



giving the cluster number density within the comoving volume 
element dV/d zdQ for a given solid ang le AQ on the sky. For the 
mass function lPress & Schechted (Il974i) proposed 



dn 
dM 



(M,z) = 



2 po^c d In ctr 
no-RD+(z)M^ dlnM 



exp 



2crlDl(z)j 



(11) 



which gives the number density «, of halos in [M, M + dM] at a 
given redshift z, where po is the present matter density, cr^ is the 
variance of the mass fluctuations and D+(z) is the growth factor. 
We utilise the Press-Schechter mass function because it does not 
depend on numerical simulations that assume a certain cosmo- 
logical model. By integrating over the mass and redshift interval 
of interest, one obtains the number of clusters in the aforemen- 
tioned bins on the full-sky 



N, 



AM,Az 



: 47r 



Jaz Ja. 



dN 

dM (M,z). 

AM dMdz 



(12) 



The number of clusters we assign to a bin is then obtained 
from a Poisson distribution with an average value of Nam.Az- 
In order to create the sample of the A^ clusters in the given bin 
(AM, Az) we use a uniform random distribution. The mass bins 
are equally logarithmically spaced with log (AM) = 0.1 start- 
ing at the minimum mass Mn,in = lO'^* /j^'Mq up to the upper 
limit of Mmax = 1 X 10'^ /j^'Mq and the redshifts range from 
Zmin - 0.01 up to Zniax - 3, whcrc the size of the redshift bin is 
Az - 0.02. Having obtained the mass and redshift distributions 
of the clusters, the remaining task is to model the individual SZ 
signal and the spatial distribution on the full-sky. 

4.3. Modelling the SZ signal of the Individual clusters 

After the random realisation of A^ clusters in AM and Az, both, 
mass and redshift, for the individual clusters are known. The next 
step illustrated by the right column in Fig.|2] is to model the indi- 
vidual Compton-y signal, which can be derived from the cluster 
scaling relations. The mass is connected to the virial radius ryir 
by M = 47r/3AvPciit'"vii-' where Ay is the mean overdensity of a 
virialised sphere and pcrit is the critical cosmic density. Solving 
for Tvii- gives 



9.5103 (n^Aviz) 



1 +z 



ilm(z) 



-1/3 



Mel 



1/3 



/r'Mpc. (13) 



Moreover, application of the virial theorem allows us to connect 
the mass to the temperature by means of 



ksT, ^ /3^\l + z) 



^mAv(z) 



1/3 



Md 



10'5Mo/i-i 



2/3 



keV, 



(14) 



with the normalisation /3j = 0.75, being valid under the as- 
sumption of hydrostatic equ ilibrium and isothermality (see e.g. 
lKaiseii[T986t ISaraziniri986h . However, since we want to model 
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Fig. 3. Redshift evolution of the cluster properties for the ACDM and the EDE 2, 3 and 4 models. The upper left panel shows the 
number of Poisson-distributed clusters with M > 10'^ /i^'Mq in redshift bins of Az = 0.02 out to redshift zrnax - 3. The other three 
panels show the properties of individual clusters, where solid lines denote a cluster of 10'^ /i 'Mq and the dashed lines are for a 
system with 10'"* /i^'Mq. The upper right figure shows the central Compton-y parameter yo, the lower left the total Comptonisation 
Y and lower right the virial angular extent 0vir- 



the SZ signal for a given cluster, it is necessary to model the spa- 
tial distribution of the gas component, which is given by the fi- 
profile if hydrostatic equilibrium and isothermality are assumed. 
Therefore the distribution of the thermal electron density is given 
by 



-3/2/J 



(15) 



where denotes the core radi us. The core radius an d the virial 
radius ryir are related (see e.g. lGeisbiisch et al] (l2005h ) by 



(16) 



rc(z)=^(z)rvir = 0.14(1 +Z)'/Vvir. 

Due to the lack of an according relation for EDE cosmologies 
we apply this relation also in these cases. In what follows we 
chose to set y6 = 2/3, delivering good fits to the X-ray surface 
brightness of observed clusters. 

After these prerequisites, the first step is to fix Wgo by solving the 
condition 47r ^ rif-ir) r^dr - for n^Q, which directly leads to 



"eO 



47rr; 



.3 1 



+ arctan 



(17) 



Since we do not integrate to infinity we are utilising a truncated 
y6-profile. By using the assumption of isothermality in equation 
([3j and integrating out to ryir we then obtain for the central 
Comptonisation 



yo 



- - arctan(^) 



(18) 



In addition, the line-of-sight integrated Compton-y profile is 
needed and can easily shown to be given by 



y(9) - 2 — : arctan 



1 + - 



(19) 



where 6 - rjD/^. In Fig. [3] we show the redshift dependence of 
some individual cluster properties, as well as the cluster abun- 
dance for the ACDM and EDE cases. The upper left panel 
clearly shows the expected increase in the number of clus- 
ters towards higher redshifts. The other three panels show in- 
dividual properties of the clusters, like central and integrated 
Comptonisation, as well as angular extent on the sky, where 
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10^ lo-" 
multipole order I 

Fig. 4. Angular power spectra for the spatial distribution of clus- 
ters on the sky, assuming an averaged bias scheme for all con- 
sidered cosmologies for a gaussian window function. The small 
panel in the plot shows a zoom in on the peak to illustrate the 
difference between the models. 



all of these are only mildly affected by the presence of EDE. 
Therefore the main impact of EDE is the change of the abun- 
dance of galaxy clusters. 

Now the machinery for calculating the signature of the indi- 
vidual clusters in the HEALPix map is complete. First, we cal- 
culate the positions of all pixels contained in a disc of an angular 
size of 0vir, then we assign the corresponding y(6) calculated at 
the pixel centre, and in order to assure consistency with the an- 
alytic total Comptonisation Y, we scale all pixels by a factor a 
given by 



(20) 



Here y, is y in the centre of pixel / and Qpix is the pixel solid 
angle. For clusters that are unresolved by the HEALPix tessella- 
tion, we scale Y by the ratio of cluster solid angle to pixel solid 
angle. The whole procedure is performed for all clusters in the 
given mass and redshift range and the positions of the centre 
pixels are obtained as described in the following subsection. 

4.4. Placing the clusters on the sky 

After being able to model the SZ signal for clusters of given mass 
and redshift, the remaining task is to place them in a realistic way 
on the sky. Although it is quite popular to just randomly place the 
clusters on the sky, we decided to take the cluster correlation into 
account. The scheme for obtaining the spatial distribution on the 
sky is presented by the left column in Fig.|2] where the starting 
point for computing the angular correlation (or equivalently the 
angular power spectrum) of the clusters is the cold dark matter 
power spectrum given by 



P(k,z)^Ak"T^(k,z)Dl(z), 



where D+ (z) denotes the grow th factor, T(k) is the CDM transfer 
function (iBardeen et al.lll98^ . n is the primordial power spec- 
trum index and the amplitude A can be normalised using erg via 



o-i- r 

Jo 



k^Ak^ 
2^ 



■P(k)\WR{k)\^ 



(22) 



using the window function WR{k) in ^-space. Two popular 
choices for W«(fc) are the spherical top-hat window function 
Wsik) = 3/ikRy (sin kR - kR cos kR) and the Gaussian smooth- 
ing WR{k) - exp [-(fc7?)^/2]. In our case both choices lead to 
similar results for the angular power spectra. The next step is to 
correct for the biasing b(M) between the halos of mass M and 
the dark matter distribution 

5h = ^(M)(5dm, (23) 

where (5h is the halo- and 5dm denotes the dark matter density 
contrast. For the b iasing we utilise the relation introduced by 
iMo & Whit"3 ( Il996i) . which is given by 



b{M, z)^l + 



1 



- 1 



(24) 



introducing the mass dependence on the right hand side via the 
smoothing radius R 'ma-R. In order to simplify the treatment it is 
convenient to introduce the concept of an effective bias defined 
as 



^eff(z) = 



■b{M,z)^AM 



'Mmi, 



dn(M,z) 
AM 



AM 



(25) 



which is a mass weighted integral of b{M, z) over the mass, such 
that we can write 



Pu(k,z)^blg(z)P(k,z). 



(26) 



So far, everything is defined in 3D-space, but our goal is to obtain 
the angular cluster correlation with respect to the positions on the 
sky and to neglect correlations in redshift space. In view of this, 
it is natural to work with line-of-sight integrated quantities like 



S(n) 



f 

Jo 



(27) 



which is the line-of-sight integrated density contrast in direction 
n for a window function q(x) where x is the comoving distance. 
The expansion of 6{n) into spherical harmonics gives 



6{n) = ^e'"YT(fy 



(28) 



Putting all together, and applying Limber's approximation, one 
obtains for the connection between the dark matter power spec- 
trum and the angular power spectrum the relation 



rx 
Jo 



dx' 



X" 



(29) 



(21) 



where q(x) denotes the window function, which is for our case 
a simple step function being non-zero for the intended redshift 
range. 

Having obtained the angular power spectrum it is straightfor- 
ward to create a Gaussian realisation u sing t he synalm routine 
of the HEALPix package (iGorski et al .1120051) . representing the 
probability of having a cluster in a given pixel. After normalisa- 
tion of the map, the number of clusters per pixel is drawn from a 
Poisson distribution, such that the expectation value of the clus- 
ter density equals the mean of the Poisson distribution. Figure |5] 
shows one full-sky Compton-y map for the case of EDE 4 with 
the angular correlation included. 




Fig. 5. Simulated full-sky SZ-map for the EDE 4 model, taking the angular correlation into account. The colour shading indicates the 
Compton-37 parameter, being proportional to arcsinh(10''3'), such that the maximum value of 7. 167 corresponds toy - 6.48 x 10""^. 



5. Simulating SZ observations by Planck 

In order to simulate realistic Planck observations one has to take 
several complications into account that are briefly summarised in 
the following subsections. First, the CMB is not the only emit- 
ting source in the frequency range of interest, our Galaxy for 
instance, radiates through several mechanisms like synchrotron, 
dust and free-free emission. Second, the detectors suffer from 
instrumental noise, and the scan-pattern as well as the frequency 
response in the different observing channels have to be taken into 
account. Planck will operate in nine frequency channels ranging 
from 30 GHz to 857 GHz with an angular resolution starting at 
33.4 arcmin for the lowest frequency channel and going down to 
5 arcmin for the high frequency channels above 143 GHz. The 
foreground components and the instrumental noise u sed for the 
present study are the same as in Schafer et alj (l2006al) . 



5. 1 . Foreground components 

In order to perform realistic simulations of Planck observations 
it is necessary to include several foreground com ponents, being 
active in the frequency range of interest. Following lSchafer et alJ 
(1200 6b') we assume isotropy of the spectral emission properties, 
as well as that the emission amplitude can be modelled by tem- 
plate extrapolation. For our studies we included the compulsory 
CMB, as well as Galactic dust, synchrotron and free-free emis- 
sion: 

- CMB: For both the AC DM and the EDE cosmologies 
CMBEASY (iDoranI 12001) was used to generate the Ce co- 
efficients of the CMB power spectra. From the right panel of 
Fig. [T] it is evident that the Q's for the different models are 
indistinguishable on cluster scales. The corresponding sets 
of flfm's have been generated using the synalm routine. 

- Synchrotron: Being dominant at frequencies below 100 
GHz and extending to high Galactic latitudes, the syn- 



chrotron emission is an important foregroun d component. 
The used tem plate is based on observations jHaslam et al.l 
Il98li Il982l). being adapted to Planck observations by 



Giar dino et al.l (|2002|) . The spectral slope of the synchrotron 
emission has a spectral break at v = 22 GHz, as it has been 
observed by the WMAP team. 

- Dust: For the high frequency range above 100 GHz Galactic 
dust emission is the dominant foreground c omponent, being 
mostly concentrated in the Galactic dis c jFinkbeiner et al.l 
119991 12000t ISchlegel etal]ll997l Il998f) . The template de- 
rived from observations at v = 3 THz has been extrapolated 
to the Planck channels, using a two-component model by C. 
Baccigalupi, where the emission is given by the superposi- 
tion of two Planck laws with T[ - 9.4 K and T[ - 16.2 K at 
a fixed ratio of 0.49. 

- Free-free: The template of the Galactic free-free emission is 
based on Ha-observations (Finkbeiner 2003), and the spec- 
tral model suggested by Valls-Gabaudi (11998 ) for the conver- 
sion of the Ho- intensity to the free-free intensity. The conver- 
sion can be parameterised by the plasma temperature, lead- 
ing to a v^-law for the spectral dependence of the free-free 
brightness temperature. 

By this, the list of contaminants in Planck's, frequency range is 
by far not complete. In the current analysis we omitted the con- 
tribution of Zodiacal light and the emission of Solar System bod- 
ies, as well as the contamination by microwave point sources, 
such as star-forming galaxies and active galactic nuclei (AGN) 
(Toffolatti et al. 1998). Even though the impact of point sources 
on cluster detection is for sure not negligible, we did not at- 
tempt to model point source contamination, since the spectral 
behaviour, spatial clustering and biasing, as well as the AGN 
duty cycles, are poorly known, especially for the case of EDE 
cosmologies. In this sense we assume that the microwave sky has 
been cleaned of point source contamination prior to the cluster 
survey. 
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filtering 



Fig. 6. Flowchart for the preparation of the coadded maps and 
subsequent filtering procedure. 

5.2. Instrumental issues 

Apart from the foreground contamination, several instrumental 
issues have to be taken into account when modelling Planck 
observations, such as the angular resolution (beam size), detec- 
tor noise and the frequency response of the individual channels, 
which are briefly summarised in the following: 

- Planck beams: The beam shapes for the nine channels are 
approximated by azimuthally symmetric Gaussians b{6) - 
(2no-ey^ exp(-e^/2crl) where crg^ Ag/ V81n(2) with Ag be- 
ing the angular resolution, starting from 33.4 arcmin for the 
lowest channel and going down to 5 .0 arcmin for all channels 
above 143 GHz. 

- Noise maps: The noise per pixel is obtained by drawing a 
Gaussian distributed random number with zero mean and 
and width ctn in units of antenna temperature, which is ob- 
tained from the Planck detector database. Assuming Poisson 
statistics, the noise level can be scaled down by VMiet with 
A^det being the number of redundant detectors for each chan- 
nel, and by y/Nhh, which is the number of observations in the 
pixel of interest. 

- Frequency response: Planck's frequency response can be 
well approximated by a top-hat, centred on the fiducial cen- 
tre frequency vq, with a channel dependent window size Av. 

5.3. Preparation of coadded maps 

The preparation of the observed coadded sky maps, including 
all astrophysical components and instrumental properties, is pre- 
sented in the upper part of the flowchart shown in Fig. |6l First, 
the skymixer module, being part of the Planck LevelS pipeline 
dReinecke et al. 2006), is used to add up the templates of the 
CMB, the galactic foreground components and the SZ effect in 
Ofm-space, applying the spectral laws for each component for the 



given frequency response in the respective channel. After that, 
the maps are convolved with the beam profiles of the respec- 
tive channel and the noise maps are added, leaving nine coadded 
maps C(m in flf,„-space that can directly be fed into the filtering 
routines as described in the next section. 

6. The filtering 

In order to study the impact of EDE on the observed clus- 
ter sample we decided to use the multi-fre quency-matched- 
filter (MFMF) method intro duced by ( Schafer et alj l2006bl: 
ISchafer & BartelmannI |2007|) . It is a generalisa t ion o f the 
matched filtering, as introduced by ISanz et all (1200 1[) and 
iHerranz et al.l ( |2002[) . to spherical datasets formulated in spheri- 
cal harmonic space. The derivation of the filter kernels ^{\0\,R), 
being constructed for the detection of objects of a certain scale R, 
can be formulated as the solution of a variational problem. The 
detection map R) for channel / is obtained by convolving 
the corresponding filter kernel ^i{\9\,R) with the coadded map 
Ci(0), where the variance crj^iR) of the detection map is con- 
strained to be minimal. Additionally the filter kernels have to 
fulfil two further conditions: 

1. to ensure the existence of a scale R^), such that the filtered 
field (£>(0source,^o)) at the source position ^source is maximal 
and that 

2. the filter kernel should be an unbiased estimator of the un- 
derlying source amplitude. 

The source profiles needed for the filter construction are as- 
sumed to be spherically symmetric y6-profiles and to be superim- 
posed on a homogeneous and isotropic fluctuating background, 
obeying Gaussian statistics. In this sense it is further assumed, 
that all foregrounds have isotropic spectral properties. Further 
detail on the filter construction, as well as the results of extensive 
testing of the filter perfo rmance and the ke y ass umptions m ade, 
can b e found in dSchafer et alj l2006bl; ISchafer & BartelmannI 
lIOOTh . 

For our case the filtering procedure is summarised in the 
flowchart shown in Fig. |6] After the construction of the coad- 
ded Compton-y maps in flf,„-space, as shown in the upper part 
of the flowchart, filters optimised for finding beam convolved 
;0-profiles are constructed and normalised by the variance of the 
coadded sky maps. In doing so, the convolution of the filters with 
the coadded maps and back-transforming to real-space, yields a 
detection map, giving a signal-to-noise ratio (S/N) for each pixel. 
Applying a S/N threshold of 3 o" and a Galactic cut of 10 deg, 
prepares the detection map for the last step of cluster identifi- 
cation. In this step all peaks above the threshold are identified 
and the position in Galactic coordinates, as well as the peak S/N 
are dumped to the cluster catalogue, which can afterwards be 
cross-checked against the input cluster catalogue. The fact that 
the method works on the full-sky, is fast and blind makes it an 
ideal tool for our purposes in allowing fast parameter studies. 
For a given set of coadded maps it take less than 45 minutes on 
a normal PC to obtain a cluster catalogue. 

7. Results 

To perform our study we created four full-sky Compton-y 
maps; three for the EDE 2, 3 and 4 models and one for the 
fiducial ACDM cosmology. The maps contain a cluster sample 
with masses ranging between 10'^ /z^'M© and a few times 
lO''' /z^'Mq out to redshifts of Zmax = 3.0. In order to speed 
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Fig. 7. Number of confirmed detections in bins of S/N of width A^g^ — 0.5 for Ftieshoid = 2 x 10 arcmin^ and an association 
radius of rsearch = 30 arcmin for the EDE 2, EDE 3 and EDE 4 models from left to right. The upper row shows the clean datasets 
(CTN) and the lower one the foreground contaminated datasets (CTNSDF). For comparison also the ACDM case is also shown in 
each plot. The grey strip on the left indicates the filtering detection threshold of 3 cr. 



Table 2. Properties of the simulated maps, Neon is number of 
correlated clusters and A^jot is the total number of simulated sys- 
tems. 



map 


^corr 


A'tot 


ACDM 


1.02 X 10^ 


3.58 X 10' 


EDE 2 


3.08 X 10^ 


5.68 X 10' 


EDE 3 


2.46 X 10^ 


5.02 X 10' 


EDE 4 


4.10 X 10^ 


6.21 X 10' 



up the computation, the angular cluster correlations have been 
taken into account for systems with > 1.5 x lO''* /i^'Mq, 
such that we assume a homogeneous and isotropic SZ back- 
ground. This threshold has been chosen in order to ensure 
numerical feasibility, as well as minimising the effect on the 
detected cluster sample. Varying the threshold actually showed, 
that the statistical impact on the detected cluster sample for 
our filter method is negligible. Depending on the model, the 
maps comprise in total more than sixty million galaxy clusters 
as summarised in Table |2] ensuring a sufficient treatment of the 
SZ-background due to unresolved clusters. We fed these maps 
into our observation and filtering pipelines for the clean case, 
where we just added the CMB and instrumental noise to the 
maps, dubbed hereafter the CTN case, and for the CTNSDF 
case, where we added also the Galactic foregrounds. 
Assuming a survey threshold of Fti-eshoid = 2 x 10""^ arcmin^, 
corresponding to a limiting mass M„„„ - 3.5 x 10'^ /i^'M© 



at z = 1.0, and an association radius of rsearch = 30 arcmin 
we cross-check against the input cluster catalogues of our 
Compton-y maps and obtain in this way cluster catalogues 
that are free of false detections. In this sense we assume a 
Planck cluster catalogue comprising only confirmed detections, 
where the validity and redshift of each detection stem from 
follow-up observations. Moreover, our study does not aim to 
recover individual cluster properties, like mass or the integrated 
Comptonisation, but to exploit the effect on the detected number 
counts. 



The first step is to look at the distribution of the number 
of detections in bins of S/N, as shown in Fig. [T] which is the 
information directly obtained from the filtering process. For 
the clean CTN case (upper row) as well as for the CTNSDF 
case (lower row), the curves in Fig. [T] look as expected, in the 
sense that the number of detections is a decreasing function of 
the S/N. However, the deviations between the ACDM case and 
the EDE cases are evident. For S/N above 5cr the number of 
detections is much higher (up to 100% and more) in some bins, 
and overall the distribution falls of flatter with the S/N. This 
result is not surprising since we expect in EDE cosmologies 
a higher abundance of massive systems in the past and these 
systems are quite likely giving rise to higher S/N. However, 
when having a closer look on the EDE 3 CTNSDF case in Fig.|2l 
it is obvious that for this case the difference is less pronounced 
and can hardly be distinguished from the ACDM case. 
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Table 3. Properties of the detected cluster sample for Ftreshoid - 
2x10"^ arcmin^ and an association radius of rjearch = 30 arcmin, 
^det number of detections clusters and A^true ^h^ '^'ue ones. 



map 


^(lef 


A'true 


contamination [%] 


ACDM CTN 


4221 


2914 


30.96 


EDE 2 CTN 


5272 


5117 


2.94 


EDE 3 CTN 


4796 


4495 


6.28 


EDE 4 CTN 


5651 


5560 


1.61 


ACDM CTNSDF 


3389 


2123 


37.36 


EDE 2 CTNSDF 


4124 


3957 


4.05 


EDE 3 CTNSDF 


2142 


1971 


8.41 


EDE 4 CTNSDF 


3680 


3568 


3.04 



After the assumed redshift determination by follow-up obser- 
vations, it is also possible to study the distribution of detected 
systems in redshift bins. The results for a redshift bin size of 
Az = 0.05 are shown in Fig. [8] again for all models, and the 
CTN as well as the CTNSDF datasets. In the redshift range 
< z < 0.6 the increase in number of detections for EDE 
relative to ACDM is quite pronounced for the clean case and 
also for the EDE 2 and 4 CTNSDF datasets, whereas the EDE 
3 CTNSDF case is indistinguishable from a ACDM model in 
terms of the redshift distribution. As expected the strongest 
detectability is found for the limiting EDE 4 case, resulting in 
a boosted number of detections up to z ~ 1 . When taking the 
values of OoE.sf from Table [T] into consideration, one can argue 
that with the used detection procedure EDE models having 
around 4% of OoE.sf should be detectable. 

In addition to the number counts, EDE has an effect on 
the contamination of the detected cluster sample. We define 
the contamination as the ratio of the number of true detected 
^true clusters over the total number of claimed detections A^jjej-. 
The values for all four setups can be found in Table [3] from 
which it is evident that all EDE catalogues are significantly 
less contaminated than the ACDM case. Such an effect can 
be explained by the impact of the relative noise levels for the 
ACDM and EDE cases, such that for the EDE case one has 
more significantly detectable clusters relative to the uniform 
background than for ACDM. This effect can also be mimicked 
by downscaling the noise as shown in the appendix. In this 
sense EDE helps its own detectability not only by increasing the 
number counts, but also by improving the purity of the detected 
cluster sample. 

In the introduction it was mentioned that EDE also might 
offer an explanation for the excess in angular power at high 
multipoles, as discussed by ISadeh et aP (l2007h . To show that 
this is also the case for our choice of models we calculated 
the SZ power spectra followi n g the approach as discussed. 



for example, bylC oorav (2 0001) . iKomatsu & SeljaS (l2002h and 

20] 



Refregier & Tevssier (.2002 !). In order to be comparable to 



Sadeh et al .1 (I2007h we show in Fig|9] the SZ power spectra at 
an observing frequency of v = 3 1 GHz and the redshift distri- 
bution of Cf . As expected for the power spectra, all EDE mod- 
els lie above our fiducial ACDM cosmology and for the EDE4 
model the boost in power is strongest, despite having the lowest 
CTg = 0.655, showing that EDE could in principle offer an expla- 
nation to excess power at high multipoles. The redshift distribu- 
tion of the C[ clearly shows an enhancement in the high redshift 
tail, since EDE contributes strongest at higher redshifts. Even 
though the current observational status does not allow the mod- 






^DE3 


1=20000 


EDE3 




EDE4 




ACDM 


l=300o''' \\' 









Fig. 9. SZ angular power spectra at v = 3 1 GHz (left panel) and 
redshift distribution of C( (right panel), plotted as d In Cf /d Inz 
for € = 3000, 20000, for the EDE2, EDE3 and EDE4 as well as 
the fiducial ACDM model. 



els to be distinguished, and the nature of the excess is still un- 
clear, future high quality measurements of the SZ power spectra 
will be a powerful tool for confirming, or ruling out, deviations 
from the ACDM model. 



8. Summary and conclusions 

The aim of this paper is to present the results of a study on the 
impact of the altered structure formation history in EDE cos- 
mologies on Planck's, SZ cluster sample, using the spherical 
MFMF filtering approach. By computing the SZ power spectra, 
we show that EDE indeed leads to a boost in the level of power, 
which could in principle explain the excess of power on small 
angular scales of the observed CMB power spectrum. 
In order to examine the impact of EDE on Planck's, SZ clus- 
ter sample, we introduced a method for constructing full-sky 
Compton-y maps for EDE cosmologies, taking angular cluster 
correlation into account. To do so we utilise a line-of-sight in- 
tegrated approach, using Limber's approximation, to convert the 
three dimensional power spectrum to the angular one, such that 
we neglect the coiTelation in redshift space. A Gaussian real- 
isation of the angular power spectrum, normalised to the total 
number of expected clusters, leads to the expectation value of 
the number of clusters in a given pixel. We obtain then the spa- 
tial distribution of the cluster sample by drawing the number of 
clusters in a given pixel from a Poisson distribution, such that 
the mean of the distribution equals the cluster expectation value. 
Based on the full-sky SZ maps, we simulate Planck observations 
in all nine frequency channels, including the optical setup, detec- 
tor noise, frequency response, as well as Galactic foregrounds. 
The resulting detected cluster samples have been compared to 
the fiducial ACDM case after cross-checking with the respective 
input catalogues, leading to the following results: 

1. We showed that the increased presence of galaxy clusters 
in EDE cosmologies is also reflected in the detected cluster 
sample, using the MFMF filtering technique. 

2. The possibility of detection depends on the strength of the 
EDE. For the realistic CTNSDF case of EDE 3 there is no de- 
tectable effect on the redshift distribution relative to ACDM. 
It seems that a contribution of OoE.sf ~ 0.04 is required for 
significant detection with the considered filtering method. 

3. The presence of EDE also influences the contamination of 
the detected cluster sample via the relative noise contribu- 
tion. For the EDE case the contamination is always lower 
than for the ACDM case. The effect of the relative noise level 
on the contamination is briefly discussed in the Appendix. 
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Fig. 8. Number of detections in redshift bins of width Az = 0.05 for FtreshoU = 2 x 10" arcmin and an association radius of 
'"search = 30 arcmin comprising the EDE 2, EDE 3 and EDE 4 models from left to right for the clean case (CTN) in the upper row 
and for the galactic foreground dataset (CTNSDF) in the lower one. To ease comparison the ACDM case is also shown in each plot 
as well. 



Under the assumption that Planck will be able to deliver all 
necessary information on foregrounds, it should be possible to 
detect significant deviations from the ACDM expectations, pro- 
vided that the selection functions for the applied filtering meth- 
ods are sufficiently known. And since our filtering method is op- 
timised globally for full-sky application, it can be expected that 
methods working on patches outside the Galactic plane will find 
even more significant differences than our analysis. 
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Fig.A.l. Contamination (grey bars) and completeness (green 
bars) for the ACDM CTN dataset as function of the scaled CMB 
for a Fthrestioid = 3 X 10""* arcmin^ and rsearch =15 arcmin. The 
solid lines show the according contamination levels for the CTN 
case of the models EDE 2, 3 and 4. 



Appendix A: Impact of noise levels to the filter 
performance 

In order to study the surprising results in terms of the contam- 
ination of the detected cluster sample, we decided to perform a 
toy model test by downscaling the level of the CMB, which acts 
as source of noise to the filtering process. We applied in this case 
more strict association criteria of Fthrestioid = 3 x 10"'* arcmin^ 
and Tsearch = 15 arcmin and included also the completeness, de- 
fined as the ratio of A^true o^^r the number of cluster in the map 
above Fthreshoid, in our analysis (denoted by the green bars in Fig. 
lA.ll By downscaling the CMB, one can mimic the effect of the 
relative noise level on the properties of detected cluster sample. 
The results of this test are shown in Fig. lA.U where it can be seen 
that the contamination is decreasing with the decreasing noise 
level. The coiTesponding contaminations for the EDE 2, 3 and 4 
models are illustrated by the solid lines and it can be seen that 
the levels of the EDE cosmologies can be reached by a relative 
lowering of the noise level. Hence, the presence of EDE does 
not only decrease the number of detectable clusters, but is also 
lowers the contamination by increasing the relative noise ratio 
between the background and the enhanced low and intermediate 
redshift cluster sample. 



